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A model of scalar turbulent advection in compressible flow 
is analytically investigated. It is shown that, depending on 
the dimensionality d of space and the degree of compressibil¬ 
ity of the smooth advecting velocity held, the cascade of the 
scalar is direct or inverse. If d > 4, the cascade is always 
direct. For small enough degree of compressibility, the cas¬ 
cade is direct again. Otherwise it is inverse, i.e. very large 
scales are excited. The dynamical hint for the direction of 
the cascade is the sign of the Lyapunov exponent for parti¬ 
cles separation. Positive Lyapunov exponents are associated 
to direct cascade and Gaussianity at small scales. Negative 
Lyapunov exponents lead to inverse cascade, Gaussianity at 
large scales and strong intermittency at small scales. 

PAGS numbers 47.10.-bg, 47.27.-i, 05.40.-bj 

The keystone of the celebrated 1941 Kolmogorov- 
Obukhov |Q,|| theory for 3D fully developed turbulence 
is the direct (downscales) energy cascade. Many other 
examples of direct cascades have later been found for 
turbulent systems (see for a review). The pres¬ 
ence of a direct cascade expresses the fact the average 
flux of an integral of motion, which holds for the system 
in the absence of forcing and dissipation (e.g. energy 
for 3D Navier-Stokes turbulence and many examples of 
wave turbulence, vorticity for 2D Navier-Stokes turbu¬ 
lence, etc.),is directed toward small scales and is con¬ 
stant along the scales. The variety of turbulent systems 
is not exhausted by direct cascades. When the pumping 
is supplied, the flux of an integral of motion can go to¬ 
ward large, and not small, scales. This is the case of 2D 
turbulence, where the dynamical constraints due to the 
presence of two integrals of motion lead to the remarkable 
phenomenon of the inverse energy cascade discovered by 
R.H. Kraichnan in ||]. Other examples of inverse cascade 
are known in wave turbulence ^. 

Advection of a passive scalar 9{t;r) (it might be the 
concentration of a pollutant or temperature) by a tur¬ 
bulent incompressible flow belongs to the class of direct 
cascades |^. The direct cascade survives if the Navier- 
Stokes turbulent velocity is replaced by a synthetic field 
with prescribed statistical properties (this model, 
introduced by Kraichnan in |^ , is attracting a great deal 
of attention for the anomalous scaling discovered there 
P, pR|JTl| ]). A one dimensional compressible generaliza¬ 
tion of the Kraichnan model was recently introduced in 


|]l^ . We have considered the smooth limit of the model 
in and shown that an inverse cascade takes place. 
This has led us to investigate the general relation be¬ 
tween compressibility and the direction of the cascade. 
The aim of this letter is to present and analyze a model 
where we can continuously move from inverse to direct 
cascade by varying two parameters: the dimensionality 
of space d and the degree of compressibility. 

The dynamical equation is 

(5t-f Ua(t; r)VQ-KA) 0(t; = r), (1) 

for the Lagrangian tracer scalar held 0(t; r) (say tem¬ 
perature or entropy), steadily supplied by the random 
Gaussian source (j) and advected by the compressible 
d-dimensional Gaussian random velocity u{t',r), hav¬ 
ing zero average. The correlation of the pumping is 
= x(ri - r 2 )^(ti - ^ 2 ), with x(r) reg¬ 
ular at the origin and decaying fast at distances larger 
than the integral scale L. The molecular diffusivity k 
is supposed small enough to be in the fully turbulent 
regime, i.e. L is much larger than the dissipative scale. 
The pair correlation function of the velocity fluctuations, 
{ 6 ua{t;r) 6 ui 3 {t']r)), is 

2(dC^-5'^)r„r^-(2C*^ - (d-H)5'^)daffr% ^^ 

d{d-l){d + 2) ^ ^ ^ 

where = ((VcWq)^), = {{VaUp)'^) and Sua{t;r) = 

'^a(t, r) llciit, 0). 

We show for the model that: the cascade of the scalar 
is inverse if d < A and the degree of compressibility 
C'^ IS'^ > d/A ; otherwise it is direct. The dimension d = A 
turns then out to be critical for the direction of the cas¬ 
cade. When the cascade is inverse and no infrared cutoffs 
are present, the energy is clearly piled up at the large 
scales and moments of the scalar field diverge linearly 
with time. Gorrelations of scalar differences will however 
reach stationary values. The PDF at the steady state 
of scalar differences is found explicitly in the convective 
interval, i.e. at scales much larger than the diffusive one, 
both for direct and inverse cascades. The Gaussianity 
of the scalar distribution, known to be present at small 
scales in the incompressible case ||^, emerges at large 
scales when the cascade is inverse. Small scales are then 
shown to be strongly intermittent. 
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Before the systematic analysis, it is worth presenting 
first some simple intuitive arguments. Since the scalar is 
a tracer, its statistics is very directly related to the one 
of the Lagrangian trajectories separation and, specifically 
for smooth velocities, of Lyapunov exponents, describing 
the rate of the exponential-in-time stretching (or con¬ 
traction) of the separation. We shall derive below the 
following relation between the sign of the maximum Lya¬ 
punov exponent and the direction of the cascade: A is 
positive in the case of direct cascade and negative in the 
case of inverse cascade. Its physical origin is simple to 
explain. Indeed, once the statistics of the scalar field is 
expressed in terms of Lagrangian paths properties, the 
crucial question becomes: given two particles initially 
separated by r, what is the probability that their dis¬ 
tance will ever reach the typical scale of the pumping 
L ? If the Lyapunov exponent is positive, separations 
will typically grow exponentially in time. Starting from 
a separation r L, the scale of the pumping is almost 
certainly reached ; viceversa, if r ^ L, there is practically 
no chance to reach L. This is the dynamical hint of the 
direct cascade. Consider now a negative Lyapunov expo¬ 
nent. The picture is completely reversed. Since typical 
trajectories are contracting, the probability to reach L is 
much higher at scales r ^ L than r L. The character¬ 
istic sign of the inverse cascade is emerging here. Further¬ 
more, the described physics of the relation between the 
stretching-contraction interplay does have a certain gen¬ 
erality. For incompressible flow, the strain tensor being 
traceless, no strong trapping phenomena are possible and 
the rate of typical separations growth is expected to be 
positive. Strong trapping appear for compressible flow, 
where it can lead to a substantial slow-down of transport 
(see 1^) and, possibly, to a negative rate of the La¬ 
grangian separations stretching. Finally, the importance 
of traps for transport should definitely reduce when the 
dimensionality of space increases. 

Let us now start the systematic analysis of the model 
by considering Lagrangian separations statistics. The 
stochastic differential equation governing the evolution of 
the separation Ra{t) between two Lagrangian particles, 
in the absence of molecular diffusion, is dtRa = o'apRp- 
The statistics of the random strain (t“^ is Gaussian. 
Its irreducible correlation function {{a.ya(t)asf 3 {t'))) is 
simply obtained operating with Vf.Vl/2 on (^). The 
average {uap) = is fixed by the condition 

{Ra{t)) = Ra{0), which is an immediate consequence 
of isotropy (see for more details). It is worth noting 
that the average of a vanishes for incompressible flow. 
The 2n-th moment of R{t) obeys the following differen¬ 
tial equation 


dt{R^n=^{R^n 

d 



n — 1\ 

r 

d -\- 2 ) 


2(n-l) 
d + 2 



( 3 ) 


corresponding to the following Gaussian statistics of the 
exponential stretching rate A = \n[R{t)/R{0)]/1: 
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The largest Lyapunov exponent (generally there are d of 
them) has been denoted by A and ( is the variance of A. 
For incompressible velocity fields, = 0 and the Lya¬ 
punov exponent A is always positive. The opposite limit 
is the one of gradient-type velocity fields u = where 
the equality holds. The interest¬ 

ing conclusion arising from (Q) is that d = 4 is actually 
a critical dimension for gradient-type fields. For generic 
smooth flow, the largest Lyapunov exponent is always 
positive for d > 4 and its precise behavior (including the 
value of the possible critical dimension dc = /S^) 

depends on the specific value of the degree of compress¬ 
ibility C'^fS'^. It follows from (^ that, if the Lyapunov 
exponent is negative, the low-order moments of the La¬ 
grangian separation decay in time, but their high-order 
(e.g. integer moments) grow exponentially. As first 
highlighted in , this means that the dynamics of La¬ 
grangian separations is dominated by rare events and this 
is the origin of the strong intermittency of the scalar at 
small scales which will be evidenced below. 

We proceed now with the derivation of the scalar- 
statistics. The relation between Lagrangian trajecto¬ 
ries and the passive scalar field in (^ is very direct. 
The solution of the equation can indeed be presented as 

t 

9[t\r) = f dt'(j){t'■, p{t')), where the Lagrangian trajec- 

— OO 

tory p{t) satisfies p{t') = u(t';p(t')) and p{t) = r . For 
the sake of simplicity of the presentation, we have not 
taken into account molecular diffusion (see, for example, 
Ip^ for the detailed path integral formulation including 
diffusion). Its smearing effects on Lagrangian trajecto¬ 
ries can be neglected in the analysis presented in the 
sequel. The simultaneous 2n-th order scalar correlation 
function can be rewritten in terms of the average over 2 n 
Lagrangian trajectories as 


(6*1 •• • 02n) 



permut. ^=1 o 


^ikjk (^) 

L 


( 5 ) 


where the sum is over all the possible permutations 
{ii, ...,in, ji, ■■■, jn} of the indices {l,...,2n}. The La¬ 
grangian separation I^ (t) satisfies the Langevin equa¬ 
tion associated with (® and the two Lagrangian parti¬ 
cles are initially located at and r^, respectively. T 
stands for the time of evolution and should tend to in¬ 
finity for the system to attain the stationary state. The 
statistics of the whole set of d Lyapunov exponents is re¬ 
quired to reconstruct the Lagrangian dynamics of a gen¬ 
eral d dimensional structure. However, to analyze the 
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structure functions of the scalar (or multi-points corre¬ 
lation functions where all the points are elongated on 
a straight line), we do not need to go into the heavy 
details of subleading Lyapunov hierarchy. The point is 
that, modulo the dissipation that will be discussed later 
on, the collinear geometry is preserved by the dynamics 
00- This issue actually deserves some more detailed 
explanation. For the incompressible case, it has been 
found that the transition from collinear to slightly off- 
collinear configurations is not smooth 0- This implies, 
in particular, that the subleading Lyapunov exponents 
statistics is needed to describe the correct rate of ex¬ 
ponential decay for the far tails of the scalar PDF [0| . 
Being aware of this important effect, in the present let¬ 
ter we focus on the quantities that can be analyzed by 
collinear considerations. The study of the off-collinear 
effects is reserved for future publications. Here, we shall 
therefore consider only collinear geometry = nxi, with 
all the distances large enough for dissipative effects to 
be negligible. The 2n-order structure function S 2 n{x) 
can be easily reconstructed from the 2n-th order correla¬ 
tion function (w(T, ri) • • ■uj{T,Y 2 n)) of the scalar gradi¬ 
ent = dxiO{t]Yi), in the collinear geometry. This 

is simply done integrating all the coordinates Xi from 0 
to X and the final result for S2n = {{S{x) — 0(0))^”) is 


(2n- 1)!!2" 


dt [x(e^Vd/L) - x(e^‘r/L)] 


( 6 ) 


where the averaging with respect to A is fixed by (jj). 
The regularity of the pumping fhe origin has 

been exploited in and the integrations over the xds 
have been cut from below by the diffusive scale = ^/k. 
We will indeed see that in those cases when the cutoff is 
important, the resulting dependence on is logarithmic, 
thus justifying the used regularization. The generating 
function Z{r/L,q) (Fourier transform of the respective 
PDF) of scalar differences is simply restored from (|) in 
the form of a matrix element over an auxiliary quantum 
mechanics (see for more details) with the Hamilto- 
[x(e''/Pe) — x(e’'r/L)]. Indeed, 


nian H = — 
,-T\^/{2C 


z = 




7) = 0 


, where the wave function dt 


dx -l- H 


d* = 0, with 


satisfies the Schrodinger equation 

the initial condition dt(0;ry) = . The asymmetry 

between negative and positive A clearly emerges in the 
different asymptotic behaviors of dt. For negative A, the 
nonvanishing aX rj —> —oo initial condition expjAry/^C^] 
will survive in the stationary (T —> oo) limit, since the 
potential part of the Hamiltonian ( ~ g^) is also vanish¬ 
ing at ?7 ^ — oo. Analogous consideration of the positive 
A case, fixes the same boundary condition for dt, but in 
the opposite limit rj —> -|-oo (worth mentioning that a fi¬ 
nite diffusivity is needed for this). We thus arrive at the 
Fokker-Planck equation 


2A 


2g2 


y '^y+i 1+-^j y^v--^ [x(o)-x(y)] 


Z=0, (7) 


with the following boundary conditions on Z(y, q) 


at A < 0 Z 


at A > 0 Z 


Zy 


x/e 


y^O 


*oo 


>00 


Zy 


x/e 


0 , ( 8 ) 


0. (9) 




Here, ^ should be considered as a ^-independent param¬ 
eter to be replaced by r^/r after integration of (!)• 

Note that the Fokker-Planck equation ( 0 ) is the same 
as the one found directly from taking two repli¬ 
cas, averaging with the appropriate weight and discard¬ 
ing the contribution of the molecular diffusivity term. 
If the latter term is kept, the so-called anomaly term 
appears (following the field-theoretical terminology in¬ 
troduced by A. Polyakov in [|^) and the equation is not 
closed. In turbulent systems, it is generally expected that 
the anomaly term, accounting for the effects of molecular 
diffusion onto the convective interval range (for a discus¬ 
sion of the Polyakov theory in the context of the present 
passive scalar model, see |^), does not vanish in the 
zero diffusivity limit. The absence of the anomaly term 
is therefore neither a trivial nor a general fact. Deriving 
(^ dynamically within a multi-points approach with an 
accurate fusion of points, we thus proved the absence of 
anomaly in the present model. 

Let us finally show how to get from 0 and (^ or ®) 
the scalar differences PDF, with a special accent on the 
issue of universality. The PDF is indeed not globally uni¬ 
versal, in the sense that it will generally depend on the 
explicit form of the pumping function x)^^)- ^ relevant 
question to ask is then: what is universal in the PDF ? 
There are two different universality issues which may be 
discussed in this respect. A kind of universality, typical 
of the case of positive Lyapunov exponents both at large 
and small scales, suggests that the only relevant param¬ 
eter is the flux of 0^ (x(0)) pumped into the system at 
the integral scale |^. Indeed, at A > 0 the asymptotic 
solution of (^) can be found replacing x(0) — xix/L) by 
x(0). Using (|) and the fact that ^ is the smallest value 
in the problem, one gets an algebraic form for the gen¬ 
erating function Z = (r^/min{r, L})“, with the expo¬ 


nent a = \J\ + 2g2x[0]/C^ ~ ^/C^- Calculating the 

Fourier integral for the PDF in the saddle-point man¬ 
ner (the large parameter is r/r^), one gets the Gaussian 
behavior at \S9r \ Injcc*] 


V{rlL,6er) = 


exp [-(56>2A/(4 x[ 0] ln[a;*])] 




( 10 ) 


where x* = min{r, Lj/r^. For very large values \59r\ ^ 
ln[a:*], the Gaussian behavior transforms into an expo¬ 
nential tail M,pll . The detailed form of the tail requires 
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however the knowledge of the whole hierarchy of Lya¬ 
punov exponents and cannot be calculated within the 
collinear approach used in the present letter Q at d > 2. 
The unknown exponential tail does not affect the behav¬ 
ior of the structure functions of order n <C ln[a;*], that are 
equal to (2n— l)!!n!2"'ln"'[a;*]x"[0]/X*. This kind of uni¬ 
versality (and the corresponding restrictions on collinear 
considerations) can be easily extended to the case of neg¬ 
ative Lyapunov exponents for r ^ L. To this aim, we 
replace again x(0) — xi^/L) by x(0) in ( 0 ), impose the 
boundary condition (^) and perform the saddle point in¬ 
tegral, as it han been discussed above. We arrive finally 
at the expressions for the PDF and the structure func¬ 
tions (valid for r ^ L and negative A), which are identi¬ 
cal to with the replacements of min{r, L}/rd by r jL 
and A by |A|. 

For negative A, another kind of universal behavior is 
found at scall scales r L. In the vicinity of the origin, 
the PDF depends only on the second derivative of the 
pumping x^^(O)- This part of the PDF is indeed formed 
at the largest g ^ 1, where one can keep only the first 
term of the expansion of xi^jL) in the “potential” of (|^). 
We arrive then at the following expression for the PDF 
at |(50r| ^ 1 and r <C L 


fx"[oy 

V 


-VC r 


set 


x"[oy 




A/C='-l/2 


( 11 ) 


c = r[l/2 — A/C^]/(r[—A/C^]v^)- The expression ( [ll| ) 
shows that all the structure functions S'„(r) of low or¬ 
ders — 1 < n < — 2A/C^, are controlled by & and 
scale normally ~ r”. The behavior of the PDF at val¬ 
ues larger than those where (0 holds, is not univer¬ 
sal, i.e. it depends on the whole form of x(^lL)- Tor 
very large values S9r ^ 1 and at small scales r L, 
it is however possible to show, from the expression of 
(^, ||) at g ^ 1 and y <C 1, that the far tail is expo¬ 
nential, V ^ exp(—y^Cli^^rl/A). The non¬ 

universality emerges via the dimensional constant Xj that 
depends on the precise shape of the source function (for 
detailed explanations on the asymptotics, see 11 , where 
the tail was fully determined in the one dimensional case 
for a specific form of the pumping). The form of the 
tail, along with (|Tl|), shows that the structure functions 
of order larger than — 2A/C^ all scale in the same way 
S2uir) ^ ( r/L) at r L. This collapse of the 

high-order structure functions scaling exponents should 
be contrasted with the Gaussian normal scaling observed 
for the direct cascade. These different scalings are an il¬ 
lustration of the profoundly different behaviors that have 
been found for the model below and above its threshold 
of compressibility where the direction of the cascade is 
reversed. 

Illuminating 

discussions with E. Balkovsky, G. Falkovich, U. Frisch, 


K. Gawedzki, V. Lebedev, A. Polyakov, B. Shraiman, 
P.L. Sulem, and V. Yakhot are gratefully acknowledged. 
We have benefited from the stimulating atmosphere of 
the IHES, Bures-sur-Yvette, where we participated to 
the workshop on turbulence. This work was partly sup¬ 
ported by a R. H. Dicke fellowship and ONR/DARPA 
URI Grant N00014-92-J-1796 (MC), Russian Fund of 
Fundamental Researches under grand 97-02-18483 (IK), 
and the GdR “Mecanique des Fluides Geophysiques et 
Astrophysiques” (MV). 


[1] A. N. Kolmogorov, Dokl. Akad. Nauk SSSR 30, 9 (1941). 

[2] A. M. Obukhov, Dokl. Akad. Nauk SSSR 32, 22 (1941). 

[3] V. Zakharov, V. L’vov, G. Falkovich, Kolmogorov spectra 
of turbulence I, Springer, Berlin (1992). 

[4] U. Frisch, “Turbulence. The legacy of A.N. Kolmogorov.”, 
Cambridge Univ. Press, Cambridge, (1995). 

[5] R. Kraichnan, Phys. Fluids 10 , 1417 (1967). 

[6] A. M. Obukhov, Izv. Akad. Nauk SSSR, Geogr. Geo- 
phiz. 13 , 58 (1949). 

[7] G.K. Batchelor, J. Fluid Mech. 5, 113 (1959). 

[8] R. Kraichnan, Phys. Fluids 11, 945 (1968). 

[9] M. Chertkov, G. Falkovich, I. Kolokolov, V. Lebedev, Phys. 
Rev. E 52 , 4924 (1995). 

[10] K. Gawedzki, A. Kupianen, Phys. Rev. Lett. 75, 3834 
(1995). 

[11] B. Shraiman, E. Siggia, C.R. Acad. Sci. 321, Serie II, 279 
(1995). 

[12] M. Vergassola, A. Mazzino, chao-dyn\9702014. 

[13] M. Chertkov, 1. Kolokolov, M. Vergassola, ’’Inverse cascade 
and intermittency of passive scalar in Id smooth flow”, 
preprint. 

[14] M. Chertkov, G. Falkovich, 1. Kolokolov, V. Lebedev, Phys. 
Rev. E 51 , 5609 (1995). 

[15] M. Vergassola, M. Avellaneda, Physica D 106, 148 (1997). 

[16] Ya. B. Zel’dovich, S. A. Molchanov, A. A. Ruzmaikin, 
D. D. Sokolov, Sov.Phys.JETP 62, 1188 (1985). 

[17] E. Balkovsky, M. Chertkov, 1. Kolokolov, V. Lebedev, 
JETP Lett 61 , 1012 (1995). 

[18] D. Bernard, K. Gawedzki, A. Kupianen, preprint 5/97. 

[19] A. M. Polyakov, Phys. Rev. E 52, 6183 (1995). 

[20] V. Yakhot, preprint 4/97. 

[21] B. Shraiman, E. Siggia, Phys. Rev. E 49 , 2912 (1994). 


4 







